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Abstract 

We present the general solution to the Plebahski equation for an I) space that 

admits Killing vectors for an entire SU(2) of symmetries, which is therefore also 

CN | the general solution of the sDiff(2)Toda equation that allows these symmetries. 

Ch ■ 

Desiring these solutions as a bridge toward the future for yet more general 
solutions of the sDiff(2)Toda equation, we generalize the earlier work of Olivier, 
on the Atiyah-Hitchin metric, and re-formulate work of Babich and Korotkin, 
and Tod, on the Bianchi IX approach to a metric with an SU(2) of symmetries. 
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We also give careful delineations of the conformal transformations required to 
ensure that a metric of Bianchi IX type has zero Ricci tensor, so that it is a self- 
dual, vacuum solution of the complex- valued version of Einstein's equations, as 
^ ■ appropriate for the original Plebahski equation. 
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• «L Introduction 

> ■ 

We have long been interested in the Plebahski 1 formulation for an J)-space, i.e., a 4-dimensional, 
complex manifold with an anti-self-dual conformal curvature tensor 2 that also satisfies the Einstein 
vacuum field equations. Any such space is determined by a solution to the Plebahski heavenly 
equation, a constraining pde for a single function of 4 variables. By now many different, equiva- 
lent forms of that equation have been developed; however, the two most common are the original 
ones given by Plebahski: the first form, for a function u = u(p,p, q, q), and the second form, for 
a function v = v(x,y,p,q), either of which serve as potentials for the metric via their second 
derivatives: 



U ,pp^ ,qq " jP*?^ ,QP 



1 and g = 2(u yP pdpdp + u jP qdpdq + u jq pdqdp + u^qqdqdq) , 
or (1.1) 



v,xxV, yy - v xy + v, xp + v yVq = and g = 2dp(dx - v, yy dp + v jXy dq) + 2ds(dy + v yXy dp - v jXX dq) 



where partial derivatives are indicated by a subscript which begins with a comma. This approach 
already has a long history; nonetheless, there still seems to be considerable effort being made 3 to 
better understand the structure of the space of solutions, and any notion as to the behavior of 
the "general" solution is still far from being found. 

Nonetheless, to make some progress on that problem, one looks to simplify the question. A 
standard approach is to simplify the question by looking for those metrics that admit a Killing 
vector. If the (necessarily skew-symmetric) covariant derivative of that Killing vector has an 
ant i- self- dual part that vanishes, then the problem has indeed been completely resolved, 4 filling 
in some part of the solution space. In this case the constraining equation for u can be reduced to 
simply an appropriate solution of the 3-dimensional Laplace equation. 4 However, when that anti- 
self-dual part is non-zero, the Killing vector reduction instead gives a single pde for an unknown 
function O = Q(q, q, s) of only three remaining (complex) variables, which we refer to as the 
sDiff(2) Toda equation 5 because of its symmetry properties. The equation, again along with the 
form of the metric 7 that it generates is the following, where we use (p as the coordinate along the 
flows generated by the Killing vector: 

+ (e") ss = , g = V 1 + V~ l {d V + u) 2 , 

(1.2) 

V = ±fi s , 7 = ds 2 + 4e n dqdq , u = ^{Vt^dq - £l,qdq} 

This equation has been of interest in general relativity in various contexts, as well as some other 
fields of theoretical physics, for over twenty years. Many explicit solutions have been found 8 al- 
though few solutions of general type are known, especially for cases without any additional Killing 
vectors. Hoping to understand methods to find solutions without additional Killing vectors, there 
have, fairly recently, been several deliberate searches for so-called "non-invariant" solutions of the 
sDiff(2) Toda equation. 9 

However, it is true that the use of symmetries continues to be the most efficient method 
we have for solving nonlinear pde's. Fairly early in the studies of this equation, M.A. Saveliev 10 
used it as a platform to begin his study of continuum Lie algebras, 11 and even presented a 
form which, in principle, gave the desired "general solution" of the equation in terms of some 
initial conditions. Unfortunately that form is very complicated and does not seem to be useful 
for obtaining manageable and interesting, specific new solutions. Various explications of the 
symmetries, and the generalized symmetries have been made. 12 It is hoped that the detailed 
characterization of those symmetries facilitates both their use for finding additional classes of 
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general solutions as well as checking for noninvariant solutions. However, we propose a different 
approach in this article. 

There is one particular subset of all solutions that has truly received extensive study during 
the last 15 or so years; this is when the metric allows an entire SU(2) [or, equivalently, an SO(3)] 
of symmetries. 13 In this case Einstein's field equations are reduced to simply a system of ordinary 
differential equations, for functions of one remaining independent variable. At least the question 
of the solutions themselves has been completely answered for this case. On the other hand, the 
usual approach to this problem begins from an entirely different mechanism, built specifically on 
that large family of symmetries, and usually referred to as a solution of Bianchi type IX. Those 
solutions are formulated either via Schwarzian triangle functions 14 or Painleve transcendents. 15 
Because the mechanism used to obtain these solutions is quite different from one that would 
begin from Plebanski's equation it is not of as much use as we would like, from the point of view 
of using it as a starting point to work backwards and find more general classes of solutions of 
the sDiff(2) Toda equation. Indeed other researchers have also thought about this question, and 
work of Olivier 16 and also of Tod 17 has shown many of the details of relationships of some of the 
cases with solutions of the sDiff(2) Toda equation; however, the essential purpose of this paper is 
to make more explicit that sort of information about functions Q(q,q, s), constituting solutions 
of our equation, with the hope that this will help push forward that search for a much better 
understanding of the solution space. In particular we will begin with the most general form for 
a vacuum, Bianchi IX solution, in the form given by Babich and Korotkin, 18 who express their 
solutions in terms of elliptic theta functions. Then we will show how one may determine the 
corresponding variables for the sDiff(2) Toda equation. It is perhaps also worth noting here that 
other approaches to extending these solutions, by Iona§, 19 and by Ohyama, 20 have emphasized 
the desirability of using elliptic theta functions for problems related to this. 

II. The Bianchi IX Approach 

We begin here with the important details of the formalism usually used for investigations 
of Bianchi IX metrics, using notation that follows Babich and Korotkin. 18 The assumed three 
Killing vectors give foliations of the 4-space, as orbits of those Killing vectors, that are, at least 
topologically, spheres, so that a very reasonable approach is to use the Pauli 1-forms, {gi}f, as a 
Maurer-Cartan basis, for the assumed SU(2) symmetries. In addition we use a fourth coordinate 
normal to those surfaces, that we denote by /j,. The general form of the Bianchi IX metric, with 
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additional conformal factor, which may be resolved so that it has zero Ricci tensor and also has 
a conformal tensor which is either anti-self-dual or self-dual is usually stated as follows: 

i w{ w\ w\\ (2.1) 

d(Zi =£2 A £3 , da 2 = a 3 A £1 , da 3 = a x A £ 2 , 

where the three functions {«^}f and the conformal factor F depend only on fx. The requirements 
concerning anti-self-duality or self-duality of the conformal tensor, and the vanishing of the Ricci 
scalar, are turned into (a system of ordinary) differential equations that must be satisfied by the 
three functions {wi \ i = 1, 2, 3}. The vanishing of the quantities labeled as W+i below cause the 
self-dual part of the conformal tensor to vanish, causing the resultant space to have an anti-self- 
dual conformal tensor, while the vanishing of the other set, W_i cause the anti-self-dual part to 
vanish, with the opposite conclusion, where we do notice that the one version may be changed 
into the other simply by a change of the sign of the independent variable, //, its derivative being 
indicated below simply by a prime: 

W±i = =F a±i + a±i(a±j + a± k ) - a ±j a± k , 

w'a WjW k 

2a ±i = Y ±j + Y± k - Y ±i , Y ±i = ±— + — — = a ±j + a ±k , 

Wi Wi 

where we use the word "cyclic" with a fairly standard meaning, i.e., to mean that the indices 
k} should always take distinct values and in cyclic order, so that they imply each of the three 
possible sets of values 1, 2, 3, and 2, 3, 1, and 3, 1, 2. A complete derivation of the provenance of 
these equations is given in Appendix I. 

Since we have chosen, following historical precedent with the Plebahski approach to these 
problems, to concern ourselves with the anti-self-dual solution, we will ask that the W + i should 
vanish. With that choice we will write out more explicitly the coupled set of six equations which 
must be solved, and will suppress the +- subscript on what was called a +i above since this is the 
only one with which we will be concerned: 

= W+i =^ a- + djdk = di(aj + ah) , 1 

> k = 1, 2, 3, cyclic. (2.3) 
w\ + WjW k = Wi(a,j + a k ) , J 

The general case of a solution of these equations, with no constraint on the conformal factor, 
F = F(pt), will have a non-zero, although traceless Ricci tensor. The equations for the three 
functions {ai}f =1 are often referred to as the Halphen system, 21 and there are several different 
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i, j,k = 1,2,3, cyclic, (2.2) 



known forms for the solution, involving Painleve transcendents of type III, 17 , Schwarzian triangle 
functions, 14 or complete elliptic integrals. 22 We will not often need to use the explicit solutions; 
nonetheless, we prefer the form given by Babich and Korotkin, 18 in terms of theta functions. The 
general solution for the three functions {a^ | i = 1,2,3} is a 3-parameter one given in terms 
of the general Mobius transformation of the upper-half of the complex r-plane when it acts on 
the following (generic) particular solution of the system. We replace the usual variable r for the 
theta functions, which must have positive imaginary part, by r = ifi, so that we may treat \i as 
real-valued and positive, befitting its role in the form of our metric. The particular solution of 
the system is given by: 

+00 



d 

a,i = 2 — log#5_i , z = l,2,3; where < 
dfx 



d 2 = Co 2 ](0 \ifji) = e"*" e~^ m+1 ^ , 

m= — oo 

+00 

3 = C](O|iAO= e-" mV , (2.4) 

m= — oo 
+00 

tf 4 = tf[v 2 ](oM= (-ire-™ 2 ^, 



The Mobius transformation sends r — > (ar + b)/(cr + d), with ad — be = 1. More details of the 
transformation and how it scales the functions a^, are given in Appendix III. As well we describe 
the general theta functions, which also depend on a complex variable z. When those functions 
are evaluated at z = 0, as they are above, they are often referred to as theta coefficients, and are 
related to the complete elliptic integrals, K(k) and E(k). 

The functions Wi(p) are related to these others. For the case when the Ricci scalar vanishes 
they are related in a very simple way, involving the conformal factor. We are interested only 
in the pure vacuum case, so that the entire Ricci tensor vanishes. We will explain in detail in 
Appendix I that this can be done provided one chooses 

F(fu) = c (fU + do) , and Wi = a { + ^- log F , (2.5) 

for arbitrary constant values cq and do as discussed in the paper of Babich and Korotkin, 18 with 
reference to work of Tod. 25 It is also worth noting that the first solution of this sort, determined by 
Atiyah and Hitchin 22 , and explained more fully by Olivier 16 as regards its relation to Plebanski's 
equation, corresponds to the limit where cq — > at the same that Co do — > a non-zero constant, 
corresponding to a constant value for F, and therefore Wi = ai. 



III. The Matching Process 

To create the desired mappings between the two sets of variables, we will need parametriza- 
tions of the Pauli 1-forms, using a set of Euler angles on the sphere, and identified in such a 
way that the coordinate angle cp represents the variable along our original Killing vector, the one 
required by the sDiff(2)Toda equation itself. We choose one such representation given by Tod: 17 

£1 + i%2 = e l ^(d6 - i sin 9 dip) , g 3 = dip + cos 9 dip . (3.1) 

We begin the process of comparing the two forms of the metric, and establishing coordinate 
transformations by first considering that the terms involving dip 2 are uniquely picked out, since 
it is our "obvious" Killing vector. Expanding out the metric in terms of these Euler angles, we 
have the following identification, as a first step: 

1 Waia 112 772 n W\Wz 2 W 2 W 3 . 2 . 2 W!W 2 2 /n F ^ *J 

— = \\o/o(p\\ —F\\ cos ip H sin ip) sin 9 H cos 9 \ = M, , 

V 11 1 ^" u w 2 w 1 Y! w 3 J (3.2) 

M = [^3(^2 V 7 + w \ COs2 sm2 # + w l w 2 COs2 ^] ' 

where we have defined the quantity M for convenience since it will appear in many different places 
in the discussion. As the 3-dimensional portion of our metric, 7, does not depend on <p, we may 
next identify the 1-form to in the Plebanski form of the metric: 

1 f / W2W3 WsWi\ . „ . , ,„ W1W2 „ , ,1 ,„ 

w = J F 2 J — ) sin sin ^008^0+—^ cos 0dV>? 3 - 3 

v L \ Wi W2 J W3 J 

With these forms we may now calculate y(dip + to) 2 and remove that term from the Bianchi 
IX form of the metric, providing us with the form of 7 in the those coordinates: 

7 = F 4 |Mrf^ 2 + [w% sin 2 9 + (wf sin 2 ip + w\ cos 2 ip) cos 2 9} d9 2 

+ 2{w\ — w 2 ) sin ib cos ip sin 9 cos 9 d9 dtp + (w\ cos 2 ip + w\ sin 2 ip) sin 2 9 dtp 2 ^ ( 3 -4) 
= l^d[? + j ee d9 2 + 2^ei>d9 dip + ^^dib 2 = ds 2 + 4 e" dq dq . 

Our next step is to determine the factor ds 2 , which is related to the 1-form to as follows, beginning 
from Eqs. (1.2): 



2V =O jS , to = ^{Q., q dq - Q,qdq} =>- —i dto = V q ds A dq - Vqds A dq + Q, q qdq A dq , 

(e n ) o o ( 3 - 5 ) 

=S> *cfe = ^-^ds + V q dq + V q dq = Vqdq + Vqdq + (V s + 2V 2 )ds = V 2 d(2s - 1/V) , 

7 e 



where we have used the fact that Q must satisfy its constraining pde, and we note that the Hodge 
dual in question is the one generated by the 3-metric 7 in the Plebanski form, which says 



*(dq A ds) = idq , *(ds A dq) = idq , *(dq A dq) = — §e n ds . (3.6) 



As we have w and 1/V in terms of the Euler angles, we may solve this for the desired d(2s): 

(3.7) 



d(2s) = V~ 2 *tj + d(l/V) = d(l/V) + * 

7 7 



where it is re-expressed in this second form since it then involves only those quantities that we 
already know. Of course, in order to do that we must first determine the dual mapping, but 
now relative to 7, in terms of these coordinates — rather more complicated than the other ones 
because of its off-diagonal terms. We give the results below of the duals of the three basis 2-forms, 
sufficient to determine what is wanted. This calculation is straightforward, if rather lengthy, and 
is described in more detail in Appendix II, with the following result: 

ds = F 2 |[w 3 (wi cos 2 tp + W2 sin 2 ip) sin 2 9 + W1W2 cos 2 9]d/i 

+ (1V3 — w\ sin 2 ifj — W2 cos 2 ip) sin 9 cos 9 d6 + (1V2 — w\) sin?/; cos ip sin 2 9 dip^ (3-8) 
= Fdfi + gd9 + ndi(j . 

We need to square this and subtract appropriately from 7 so as to determine the remainder of 
the transformation equations, which will be done shortly. Nonetheless, we will first note here that 
this equation can be explicitly integrated, to give this Plebanski coordinate as a function of the 
ones used in the Bianchi IX approach: 

s = s(/j,, 9, VO = \F 2 [wi + W2 + (ws — wi sin 2 ip — W2 cos 2 ip) sin 2 9] . (3.9) 

Returning now to the forms given in Eqs. (3.4), we want to have 7 — ds 2 in the form 
e Q (2dq)(2dq); i.e., we want a pair of coordinates q and g, and a function O such that this equation 
would be satisfied. However, what we actually have is the difference of those two second-rank 
tensors, which we want to describe in the form above. Therefore, it is first useful to factor that 
difference into a pair of (complex- valued) 1-forms, which we call e + and its "conjugate," e_, such 
that 

e + e_= 1 -(dsf . (3.10) 



We may then find a function e n so that the separate 1-forms that make up that product may 
each be integrated. Therefore, we define a complex-valued scalar function, / = f(fj,,6,ip), to be 
chosen so that 

K e f e f = e n , =^ < _ (3.11) 

d(e- f e_) = 0j [e_= e f 2dq. 

where q and q, are the desired new coordinates. We begin with the desired factorization into the 
product of two 1-forms: 

e + =F 2 sin(6>) { [ws(w\ cos 2 tp + w<i sin 2 ip) — W1W2] cos6> + iws(w2 — wi) sin tp cos ip } d/i 
— F 2 {(wi sin 2 ip + W2 cos 2 ip) cos 2 + W3 sin 2 9 — i(w2 — w\) sinip cos ip cos 0} dO 

(3.12) 

+ F 2 sin9 {(w2 — wi) sin ip cos ip cos 9 — i(w\ cos 2 ip + w 2 sin 2 1/>)} 
= Ad/j + T d6 + Cdip , 

while the other one, e_, is obtained from e + simply by changing all the i's above to — i's, i.e., 
treating it as if all our variables are real- valued and they should be complex conjugates of each 
other. It is a lengthy but straightforward calculation to verify that this pair does indeed accom- 
plish the desired factorization scheme. The necessary "integrating factor," i.e. the function e^, 
is not uniquely determined; nonetheless the choice that we find acceptable is given, conveniently, 
in terms of its square as follows: 

e 2f =2F 2 [(w 1 - w 2 )(simjj + i cosip cos 6) 2 + (w 3 -w 1 ) sin 2 6] , (3.13) 

with the conjugate e 2 ? again obtained simply by changing the i above to —%. It is true that this 
is the square of the desired factor, rather than the factor itself; however, because of its complex 
nature it is better to display it in this form rather than insisting on just which square root is 
appropriate. It is then again straightforward algebra to show the following: 

e 3/ d(e~ f e + ) = e 2f de+ - \ de 2f A e + = , (3.14) 

which guarantees the existence of q, and also q, so that we may now write explicitly their 
differentials. 

While the necessary forms have indeed already been written out explicitly above, in Eqs. (3.12), 
the integrations that need to be performed, to determine the explicit form of q and g, are much 
more easily performed in a slightly different set of coordinates. We will choose a new set, {//, 6,p\, 
replacing ip in favor of p = i(ip + tt/2) + logtan(0/2), which gives dp = idip + d6 / sin(0). This 



will allow us to determine all the dependence on 9 quite explicitly and easily, leaving differential 
equations to be solved only in terms of n and p. We must then re-express e + in terms of {//, 6>, p}: 



e f d(2q) = e + = Ad^i + B d9 + C dp , u = coshp ; 

A = F 2 sm(9)^w 3 (wi — w 2 )u ucos9 + \/ u 2 — lj + w 2 (w 3 — u>i) cos6> j , 
B=-F 2 sin 2 9[{w 1 -w 2 )u 2 + {w 3 -w l )] = -\e 2f = -\P 2 sm 2 9 , 
C = F 2 sin(0) (wi — w 2 )u (u + a/-u 2 — 1 cos 9 j — wi 



(3.15) 



where we have given a simple symbol for the oft-repeated quantity u = coshp, and also given a 
separate name, P 2 (//,p), to that portion of e 2 ? independent of 9. The compatibility of these pde's 
for 2q has already been shown, i.e., we know that such a q exists. The integration procedure is 
explained in Appendix III, with the following result: 

2q = ±Pcos9 + N(p,p) , Psin92dq=e + , P 2 = 2F 2 [( Wl - w 2 )u 2 + (w 3 - wt)] , (3.16) 

where N = N( / u,p) is shown, in that appendix, to have several, equivalent, forms in terms of 
elliptic integrals: 



v / 7rN(^,p) = c 



d 

\(p + do)—\og'&4 L (z I ip) + nz 



z=±F(u,k)/K(k) 



F(u k) 

c (p + d ) [K{k)E{u, k) - E{k)F{u, k)] + fc , u = coshp 



(3.17) 



IV. Passing to the sDiff(2)Toda equation 

We have determined the desired coordinate transformation, in the direction 
s = s(/x,6>,V0 , "I 

q = q[jJL, 9, p(ip)} , > and also O = 6,ij)) = \ log (e 2f e 2 ~ f ^ . (4.1) 
q = q[p,6,p(il>)] , J 

However, what we wanted was O = Q(s,q, q), but the equations given are (seriously) transcen- 
dental, so that they cannot be explicitly solved to provide O = Q(s, q, q). Instead we may look at 
this set of equations as a parametric approach to that question: all the different functions named 
in Eqs. (4.1) give the desired functions in terms of "parameters," 9, and p [or p(ip)]. (We do 
note that while the integration for q was much simpler in terms of the complex-valued variable 
p, the reality conditions for the more general process are simpler to follow if we retreat back to 
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the use of ip instead.) In order to proceed from there, to show that the function O determined 
in this way does actually satisfy the sDiff(2)Toda equation, given in Eqs. (1.2), we only need to 
determine the appropriate partial derivatives, as, for instance, in the following simple example: 



d „ f da d 89 d dib d , ^ 
— = J- 1 1 — — ) O . 

ds yds d\x ds 89 dsdi/j / 



(4.2) 



Of course we also do not have the partial derivatives in the directions given just above, such 
as However, the 1-forms that form the basis for our 3-space, in Eqs. (3.12) for dq and dq 
and Eq. (3.8) for ds, along with the form for e 2 * given in Eq. (3.13), give us the entries for the 
Jacobian matrix for the coordinate transformation between these two sets of coordinates, in the 
opposite direction: 



d(q,q,s) 



(JL 6 ip 



q 
q 

s \ T 



(4.3) 



n J 



The inverted partial derivatives we need are simply the entries of the inverse matrix to this 
one. As a first step toward determining them in some convenient way, we note that it is even 
hopeful that when the rather complicated values for all these quantities are inserted into this 
matrix, Maple finds that its determinant is rather simple: 



det(J) = -| e - / e- / sin^F 6 M = -| e - n v / det7 , 



(4.4) 



where we have used Eq. (A2.3b) to insert the determinant of the 3-metric. Then the inverse of 
the Jacobian matrix itself, determined by the computer algebra program Maple, is given here: 



J" 1 EE 



d(q,q,s) 





q 


q 


s 




( f 

e J a 




c\ 


9 






V 






J\ 


J 



1 

MF2 ' 



(4.5) 



with their values given below, and we have recalled M from Eq. (3.2): 

10 



e ^MF 2 -^ = a = sin j ^3(1^1 cos 2 ip + w 2 sin 2 ip) — wiw 2 



cos0 



M = wl(wf cos 2 ip + wj sin 2 ip) sin 2 + w\w\ cos 2 

dq~ 

— iws(w2 — wi) sin -0 cos -0 1 
e~ J MF — = r = —w 3 (w 2 sm ip + w 1 cos -0) sin 6* — ti)iw 2 (wi cos ip + w 2 sin V) cos 

+ iwiw 2 {w2 — wi) sin •0 cos -0 cos , 

dip f 

e'-^MF 2 — = A = j — (1^2 — w;i)[iuiiU2 cos 2 6* + 1^3 (i^i + 1«2 — ^3) sin 2 9] sin -0 cos ip cos 9 

+ i[wiW2(w2 cos 2 i/> + w\ sin 2 ^) cos 2 9 + w 2 {wi cos 2 i/> + W2 sin 2 ?/>) sin 2 9] j/ sin 6* , 

MF 2 ^ = C =««i«)2 cos 2 + uj 2 (wi cos 2 ip + w 2 sin 2 V>) sin 2 , 
as 

90 

MF 2 — = i] = [w 2 (wi — IU3) sin 2 ?/> + w 2 (w 2 — 103) cos 2 -0] sin 9 cos , 

(Jo 

dip 

MF 2 —- = k = (w2 — wi)[w 2 sin 2 + (wiw 3 + W2W3 — u>iw 2 ) cos 2 0] sin ip cos?/; , 

(4.6) 

where the (not-displayed) partial derivatives with respect to q are just the complex conjugates of 
those with respect to g, where all the various functions we have are considered to be real. An 
interesting result that comes from this calculation is a close relationship between those partial 
derivatives that involve fx, although we have not been able to determine any useful result from it: 

^- ,a i=^% ^ 

along with the conjugate of the second equation as well. On the other hand, with these derivatives 
in hand from the inversion of the Jacobian matrix, we may now explicitly evaluate O jS , which as 
noted in Eq. (1.2) should equal 2V. We first define a name, R, for the quantity that we actually 
have in hand, from Eq. (3.13) and also Eq. (3.15), the product of e 2 ^ and e 2 -^, 

R = e 2Cl = e 2f e 2 ~f = 4F 4 | [(wt - w 2 )(sinip + i cos^cos0) 2 + (w 3 - wi) sin 2 9] 

[(wi — w 2 )(simfj — i cos ip cos 0) 2 + (w 3 — w\) sin 2 0] | (4.8a) 

= [2F 2 sin 2 \( Wl - w 2 )u 2 - ( Wl - w 3 )\] 2 = [sin 2 | P 2 1 ] 2 . 

We may then use the chain rule for d/ds, as written out in Eq. (4.2), and verify, again via Maple, 
the (required) relationship that is given in Eq. (1.2): 
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Next we use these forms to determine a form for one of the desired second derivatives in the 
sDiff(2) Toda equation, namely the second s-derivative of e n : 

8 2 R = 2<9 s (e n d s e n ) = 2[e^d 2 e Q + (d s e n ) 2 } 
1 Lo_ (d s R) 2 ] n n ^f n d 1\ (4.9) 



2 fi 



die 



2e" 



d 2 R- 



2R 



2e n V 2 (2 

os V 



It is true that the last formulation above for our second s-derivative still has an external factor 
e n ; however, when we obtain the desired formulation — in terms of the coordinates {fi, 9,ip} — for 
the other second derivative, there will be another such factor, so that they will eventually factor 
out and not cause any difficulty. However, there is the serious difficulty that all of the partial 
derivatives involving q or q contain and e^, so that it is not immediately obvious how to extract 
them in the right format. They can in fact be extracted in the desired form, but we will have to 
consider not only the product R = (e 2 -' )(e 2 -' ) that we have been considering but also the quotient, 
S = (e 2 ^)/(e 2 -^), which will allow us to take care of factors of / and / separately: 

lo gj R = 2(/ + /), logS = 2(/-7), 

2f+J= i(31o gj R + log5) , / + 27=^(31o gj R-log5). (4 ' 10) 

We begin by keeping explicit track of the factors of that appear, writing the following: 

d q = ^ + rd ° + A£ W = ef Q > ( 4 - 10 ) 

along with its conjugate form for d/dq. There are then two apparently distinct ways in which 
the desired second derivative may be written, which, of course, must be identical since ordinary 
partial derivatives commute. The first ordering is the following: 

dgd q n = \JQ (e f QlogR) = \e J Q = \e J Q{e~ f - 21 QR) 

= K"W«-(WW + 27)} = ^{«f -i^f + i 3*f } , <4 na) 

which shows that indeed we can re-write our entire equation in such a way as to only need 
the product R, and the quotient S. The other order for the initial partial derivatives gives the 
following quantity, which appears to be different: 



d q d q Q = \d q dq\og{e 2n ) = \e f Q (e f Q\ogR) = \e f Q le f 



QR 
R 



-2 e ^ e ^)-- 2 e \ R 4 R R I s R 



(4.116) 
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Since these two expressions must actually be the same we may do two useful things with them. 
As a first check on the somewhat complicated algebra, we may first insist that their difference is 
zero; namely we must have the following equality, which is simply the statement that the partial 
derivatives themselves commute: 

i [Q ,Q ]B= f!f + (4 , 2) 

The calculation (in Maple) involves quite a large number of terms on each side; however, they 
are in fact equal, verifying that all the algebra is correct. At this point then we finally want to 
determine the desired other second derivative, which takes its most symmetric form via half the 
sum of the two expressions given above in Eqs. (4.11a-b); 

1 n j QQ + QQ)R 3 QRQR IQSQR 1QSQR\ 
W-4 e | R 2 R R + 4 S R 4 S R J ' [ 6) 

As the expression for <9 2 e n given in Eq. (4.9) also has an overall factor of e n , we may add that 
expression to this one, and divide out that overall factor, reducing the verification of the sDiff(2) 
Toda equation to the question as to whether or not that sum vanishes. It is a straightforward, 
if perhaps somewhat lengthy calculation, performed in Maple, to show that this sum does in 
fact vanish, which was the necessary and sufficient condition to guarantee that the parametric 
presentation we have obtained is in fact a solution of the sDiff(2) Toda equation. 

Conclusions 

The equations we have developed give the general solution to the Plebanski equation, O = 
Q(s, g, q), when the manifold is (locally) required to have SU(2) symmetry. Second derivatives of 
Q determine the components of the metric; therefore O does not depend on the fourth coordinate 
for the manifold, cp, since variation of it has been chosen for the direction of the explicit Killing 
vector. The solution is determined parametrically in terms of an additional set of coordinates 
for the problem, {/x, 9,ip}, so that in fact we have our three desired coordinates as functions 
of them. While the presentation of {s, q, q} as functions of these original coordinates, {fi, 6>, ip} 
are explicit in terms of 9 and ip, we are unable to invert those equations explicitly because of 
the existence within them of the theta coefficients, and the (elliptic) theta functions, which are 
transcendental functions of n, analytic for all positive values of \x. It is nonetheless useful to 
collect those equations together here, from the places where they have been derived in this text: 

s = s(p, 9, ip) = \F 2 [wi + w 2 + (w 3 — wi sin 2 ip — w 2 cos 2 ip) sin 2 9] , 

_ (5-1) 

q = q(v,9,ijj) = iPcos6> + |N , q = q(ji,9,i(>) = ±Pcos6>+ ±N , 

13 



along with definitions of the symbols involved: 



P 2 = 2F 2 [( Wl - w 2 )u 2 + (w 3 - Wl )] , P 2 = 2F 2 [(w t - w 2 )u 2 + (w 3 - w t )] , 



sin ib + i cos ib cos 9 
u = coshp = —2 ; — 



(5.2) 




The overbar is used here to indicate complex conjugation in the situation where we treat the 
variables {fx, 9, ip, <p} as real- valued. Since our goal is in fact to determine general complex-valued 
solutions, that approach is still valid in the more general case, where we note that the variable 
u may take on all values in the complex plane. The equation that determines the (potential) 
function O is then given by the following, in terms of functions of {/x, 9, ifj}: 



We have shown explicitly that this parametrically-determined function O does indeed satisfy 
the sDiff(2) Toda equation, as described in Eqs. (1.2), where the function N(/i,p) is given in 
Eqs. (3.17) while the general functions Wi(fi), that determine the dependence of the metric on 
the transverse coordinate //, are given in terms of the conformal factor F = cq(h + do) and the 
solutions to the Halphen problem, cii(fi), in Eqs. (2.5). The functions Oj(//) themselves are given 
in their most general form in Eqs. (A3. 16), which constitute a (conformal) Mobius transformation 
of the in upper complex half-plane of the simpler form given in Eqs. (2.4). 

This is therefore the most general solution of the Plebahski equation that allows SU(2) 
symmetry. It is hoped that this idea of parametric solutions to a very complicated nonlinear 
partial differential equation will be a useful one, and will engender ways to determine yet more 
general solutions, with smaller symmetry groups. 

Appendix I: Calculations for the curvature of a vacuum Bianchi IX metric 

We believe it useful to show how the requirements of anti-self-duality of the conformal tensor, 
and the vanishing of the Ricci tensor, are turned into (ordinary) differential equations that must 
be satisfied by F and the three WiS. Since the purpose of the conformal factor is to affect changes 




(5.3) 



2F 2 sin 2 9 \(wi- w 2 )u 2 - [w\ - w 3 ) I = sin 2 9 \P 
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in the Ricci tensor, it is simpler to first determine these constraints by ignoring the function F, i.e., 
to set it just equal to +1, and then afterward use it to perform a conformal transformation of the 
metric, to achieve the desired results. We also note that since we are interested in the general case 
of complex- valued manifolds, for the sDiff(2)Toda equation, any particular form for the signature 
is not essential. Therefore, we will follow a fairly standard approach, and set up a tetrad with 
Riemannian signature, when variables are considered as real-valued, for our calculations, with 
F = 1: 

= (a; 1 ) 2 + (a; 2 ) 2 + (u 3 ) 2 + (a; 4 ) 2 , 

F=1 , a- (All) 

L0 4 = Hdn, co l = H—, i = l,2,3, H = y/wfuwoi , 



g 



where we use an upper-case index along with a lower-case one when we need the same values, but 
to indicate that no sum on the values is intended. 

To consider self-duality, of 2-forms, we also need a basis set for the vector spaces of 2-forms; 
we create the following two sets of triples 

E k ± ee u l A u j ± uJ k A uj 4 , i,j,k from 1,2,3, and cyclic, (A1.2) 

where the ones with the plus sign are a basis for self-dual 2-forms, and those with the minus sign 
are anti-self-dual. We use the word "cyclic" with a fairly standard meaning, i.e., to mean that 
the indices {i, j, k} should always take distinct values and in cyclic order, so that they imply each 
of the three possible sets of values 1, 2, 3, and 2, 3, 1, and 3, 1, 2. 

The curvature, and the connection, of the manifold split into separate constituents for the 
self-dual and anti-self-dual parts: the connection for the two parts, which are determined by 
separate triplets of 1-forms, along with separate triplets of 2-forms to determine the curvature. 
The complete Cartan relations are of course written as follows: 

duj a = & a r% , = dv\ + r% a r% ee \r» vXp k a a^. (ai.z) 

However, the separation into self-dual and anti-self-dual parts splits the connection and curvature 
into a pair of triplets of forms for each of them: 

G l ± ee r> ± T l4 , Qi = Q jk ± Q*4 = dG l ± - Qi A G± , (All) 

where again the upper subscripts are for the self-dual parts and the lower ones for the anti-self- 
dual parts. Inserting our tetrad, along with considerable calculation, gives us the two triplets of 
connections in the following explicit form: 

Gi=feV, 0±3 .+a ± ^±<+^, 041.5) 

\ H J Wj Wi 
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where the new functions a±k(y) are defined by the above triplet of equations in terms of the orig- 
inal tUj's, and the prime denotes the derivative with respect to [Note that within this approach 
the distinction between self-dual and anti-self-dual may be switched simply by switching the sign 
of the coordinate \i\ An immediate observation is that several of the connection coefficients are 
zero, since each of the two triplets of connections is spanned only by one basis vector, so that we 
have only 6 connection coefficients, instead of the maximal possible number of 24. These con- 
nections are then used to generate their respective curvatures, either self-dual or anti-self-dual, 
which gives us the following: 



0± = (Rjkjk ± Rj^)u J A 0J K + (R jki4 ± Riui)^ Aco 4 = Z ±l u J A u K + JC ±1 uj' A uj 4 
1 



1 \ ( H 2 \ , k fa ±I V t A (A1.6) 

< I —a : , - a ±J a ±K J w J Aw - w l I — ^ 1 a" - 1 



Once again, since each of the Q± are spanned by only two of the basis 2-forms, namely w- 7 A y fc 
and uj 1 Aw 4 , with all other possible coefficients zero, there appear to be only 3x3 = 9 non-zero 
components of the curvature tensor, namely Rjkjk, RjkiA, and Riua- However, one of them is 
not linearly independent, since the first Bianchi identity causes those three of the form RjkiA to 
sum to zero, i.e., -R1234 + -R2341 + -R3412 = 0, so only 8 of these are independent. There are also 
constraining relations between the convenient labels Z±i and IC±i. The most important of those 
are the following: 

3 3 

z +l - z_ % = ic +l + /c_, , - = = E ( /c +* + K -*) • ( AL7 ) 

i=l i=l 

To divide these components further, we separate those 8 components into those that con- 
stitute the conformal tensor, C^ u \ v , and those that define the Ricci tensor, 7£ Mi ,, and its trace, 
TZ: 

C ' avXr) RfiisXri o ^.9 i^X^-vq 9iA,r{R-v\ "I" 9vr{TZ^,\ QuX^-^rj) « {duXdvri 9u-n9v\) TZ i 

A _ fits (AL8) 

It is simpler, and equivalent, to divide the conformal tensor components into their self-dual and 
anti-self-dual parts on the basis of their definition as 2-forms, rather than their other pair of 
indices, using our two basis sets for 2-forms to accomplish this: 

Q±i±v = C i ± ilv E l ± , C±ij = C±k4: , (A1.9) 
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where the second set of equalities is a generic statement about the symmetries of this tensor, 
which in the most general case only has 5 components for the self-dual part and another 5 for 
the anti-self-dual part. However, in our particular case, using the forms for the curvature tensor 
given in Eqs. (Al.6) we find that only the following two sets of triplets of components, Ci±jk, are 
non-zero, i.e., those where all of i, j, and k are different. However, as well the sum of each triplet 
is also zero by the first Bianchi identity, so that we have only the two independent components 
for the self-dual part and another two for the anti-self-dual part: 

2H 2 C i±jk = 2W ±l - W ±j - W ±k , 

W ±l = Z i± ±lC i± (Al.lOa) 

= H 2 [Rj k j k ± 2Rj ki4: + Ri^] = =Fa±/ + a±i(a±j + a±k) — a±ja±k ■ 

Once again only two of the three elements in either one of these two triplets is independent, since 
it is straightforward to see that 

C 1±23 + C 2± 3i+C 3± i2 =0. (A1.106) 

The Ricci tensor of course involves components from both the self-dual and the anti-self-dual 
sides, which can most conveniently be described in terms of both the W±i and also the additional 
parts /Q±: 



2H 2 TZ U = - (a +i - a-i) 1 + (a +j + a-j)(a +k + a_ fc ) 

= 2 (w +j + W +fc + 2/C +l - H 2 TZ 44 ) . 



(Al.ll) 



We intend to require that the conformal curvature be anti-self-dual, which gives us only 
two constraints on the three W+i's; in order to obtain a third one, we require as well that the 
scalar curvature, TZ, vanish, which then requires all three of the W+i's to vanish. This system 
of three equations, for the three unknown functions a+i, is usually referred to as the Halphen 
system, 21 but also by various other names, including often Darboux and/or Brioschi, since there 
was considerable interest in their solution in the late part of the 19th century. Once we have 
distinct forms for the a+i's then Eqs. (Al.5) constitute a triplet of first-order differential equations 
to determine the form of the three w^s which must also be solved. A general solution of these 6 
equations has been studied by several authors; we prefer the particular form of the solution used 
by Babich and Korotkin 18 . In the case where the Ricci scalar vanishes and all three of the a + ;'s 
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are distinct they show that the difference Wi — a+i is independent of i. Referring to that difference 
as g = g(fi) = + qo), with q an arbitrary constant, we will now use that information to 
greatly simplify the previously-given forms for the traceless part of the Ricci tensor. We can use 
the ode's that we have assumed satisfied to re- write the auxiliary quantities JC+i, in a very simple 
form: 



K / W Wj w k ,WjW k 
JC+i = -a +i + a +i — = a +j a +k - a +i = {wj - g){w k - g) - [Wi - g)- 



W T Wi 



2 ( WjWk \ , / / \ 

= 9 -giwj + wk — I = -g{g + wjwj) , 

This then allows simple expressions for the components of the Ricci tensor: 

3 3 / 3 \ 

H 2 TZ M =J2lC +l = -g J2(9 + w'Jw,) = -9 I 3<7 + fl ~ ) = ^9 + 2H'/H) , 

i=i i=i V i=i Wi ) 

H 2 K U = (W +j + W +k + 2/C+, - H 2 K AA ) 

= [JC +i + g(3g + 2H'/H)] = g{g + 2[log(^/^)] / } • 



(A1.12) 



(A1.13) 



Since the sDiff2(Toda) equation generates metrics which have zero Ricci tensor, and the 
standard form of the Bianchi IX metric described above still has non-zero diagonal components 
of that tensor, we must now implement the conformal transformation of the metric generated by 
our function F = F(/j,), to arrange for the Ricci tensor to vanish as well. We consider the re-scaled 
metric, and correspondingly re-scaled tetrad basis 1-forms, denoting the re-scaled tensors with a 
"hat" over the relevant symbols: 

^^^^g^'g^V^®^, =^ £ M = W- (ALU) 

s s 

Such a transformation will leave invariant the tetrad components of the conformal tensor, C a ^5. 
However, the general transformation of the tetrad components of the Ricci tensor and, separately 
and usefully, its trace, 1Z, is given as follows: 

Kps = Kp S + Xp S , n = F~ 2 n + X , 

Xp S = F (25^1 + g p6 9 a<: ) V^cF" 1 - 3F 2 g ps [^(V^- 1 )^^" 1 )] , (41.15) 
X = -6F- 3 gP 5 VpV 5 F . 

Since we want to maintain unchanged the current zero value for the Ricci scalar, we see that this 
is straightforward provided the function F is "harmonic." Now we show that a conformal factor 
which depends only on fj,, therefore not disturbing our symmetry, is sufficient to annul the Ricci 
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tensor as is desired. Such a dependence of course simplifies the equations greatly, giving us the 
following: 



1 



F" fF'V F'H' 
3— -3 — -2 



F \F FH 



^3 ~ TJ2 ^3 ' " + 



I I p" /F'\ 2 F' ) F" 

[ T ) +2-pog(H/«, 7 )]'| , X = -6 



X^i — , 

(41.16) 



jj2 3 \ F \F I F 1 f ' if 2 F ' 

In this case the requirement that F be harmonic simply reduces to the requirement that it 
be a linear function of fx: 

= XocF" => F = cfi + d, (41.17) 

where c and d are constants. We next go to the earlier-determined, non-zero forms for the 
Ricci tensor components themselves, from Eqs. (Al.13), to see if this form for F will allow the 
transformed (trace- free) Ricci tensor to vanish. Beginning with TZ^, where we are now including 
the requirement that F" = 0, we have 

H 2 F 2 1Z 44 = H 2 7Z 44 + H 2 X 44 = -g(3g + 2H'/H) + (F'/F) [3F'/F + 2H'/H] . (41.18) 

For this to vanish we need only to require that F'/F = g: 

F = Co(/J + d ) =}► F'/F=—^—, but g=l/(jji + q ). (41.19) 

A* + do 

Therefore, by setting the two previously arbitrary integration constants, qo and do equal to each 
other, we have accomplished what was desired for TZ 44 . Since everything is now determined, we 
must now hope that this will also allow the remainder of the transformed components of the Ricci 
tensor to vanish: 

H 2 F 2 1Zu = H 2 1Zu + H 2 Xa 

(41.20) 

= g(g + 2H'/H - 2w'Jw I ) - (F'/F) [(F'/F) + 2H'/H - 2w'Jw I ] . 

We see that, yes, once again the choice that F'/F = g, is sufficient to cause these components 
to vanish as well, which tells us that the metric g = (co/g) 2 g should also be capable of being 
generated by a solution to the Plebanski equation, where g is the Bianchi IX metric described in 
Eqs. (2.1). 

Appendix II: The (3-dimensional) Hodge dual in the spherical coordinates 
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We first recall that the Levi-Civita tensor, in 3 dimensions and with an arbitrary metric g, 
with Riemannian signature, is given by 



g = det(g ab ) ?] aia2Q3 = e[a u a 2 , a 3 ] 



(A2.1) 



where the e-symbol is the usual Levi-Civita alternating symbol that takes on only the values ±1 
and 0. Then the general form of the dual of a 1-form or a 2-form, under our metric 7, is given by 



OS = \ct ah u a Au l 



*(«) = \l ac l bd Uabr]cdfU f , 



(A2.2) 



Our metric is given explicitly in Eqs. (3.4), but is perhaps best here presented in the following 
symbolic form: 

At 6 v 

1" / 7mm ^ 
7 = ((7y)) -^6 799 7^ > (A2.3a) 

ip V 7^ 7^/ 

along with its determinant: 

det 7 = 7mm [loom - To*) = 7MM( i?4 sm 2 #7^) = ( 7w ^ 2 sin0) 2 = (MF 6 sin#) 2 , (A2.36) 

where the value of M is noted in Eq. (3.2) and in Eq. (3.4) it is noted that 7 MM = MF 4 . 
The simple form of this makes it very easy to present the inverse metric: 



7 -i = ((7 «)) _ 



7a* A* 



-2 











7V></> 


-70^ 





-70^ 


700 



1 



F 2 sin0 ' ' 



(A2.3c) 



We then begin by calculating the duals of the three basis 2-forms and also, reciprocally, the three 
basis 1-forms: 



*(d9 A dip) = £dn, 

*(d/i A d6) = I Yie^dQ + ^y^dip] /7mm > 
A dtp) = -t Yieedd + w^dtp] /7 MM , 

where the various coefficients of the 3-metric, 7, may be found in Eqs. (3.4) 
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*diJ = l~ l d0 A dip , 

*d9 = i [7^^ dtp Adfj, — 70,/, d/x A d0] , (A2.4) 
*dtp = t [700 djx A dd — ^Bipdtp A d/x] , 



It is then straightforward to write down the duals of an arbitrary 1-form a and 2-form /3: 
a = H 1 dfi + H 2 d9 + H 3 dip, (3 = J 1 d6 A dip + J 2 d/i AdO + J 3 d/i A dip 



*(a) =£~ 1 H 1 d9 A dip + £(<yeeH 3 - ^H 2 )d/i Ad9 — £(^H 2 - ^H 3 )d/i A dip , 

£ £ 
*(/3) =£J x d\i H [7^2 - leeh] dO H bv>V>^2 - 70^3] dip . 

Using Eqs. (3.2-3) for 1/V and w, as well as the expressions for the coefficients of the metric 7, 
in the expression Eq. (3.7), we obtain the expression for ds as given in Eq. (3.8) in the main text. 

Appendix III: Integration for q = q(/j,,9,p), and general Theta Functions 

We intend here to show that the solution of the three (compatible) differential equations for 
q that are described in Eqs. (3.15) is in fact that given in Eq. (3.16), i.e., 



2q = |Pcos0 + N(/i,p) , Psin9(2dq) = e + = Ad^ + Bd6 + C dp , 
P 2 = 2F 2 [{wi — w 2 ) cosh 2 p — (w\ — w 3 )] . 



(A3.1) 



We begin by inserting just the first term of the form for 2q into the equation for e + . All the 
dependence on the variables is given explicitly in this case, and we obtain: 

P sin 9 d( ±P cos 9) = - ± P 2 sin 2 9 d9 + \ sin 9 cos 9 d(P 2 ) 

dP 2 

= —\e 2 f de + F 2 sin 9 cos 9 (wi — w 2 ) coshpsinhp dp + \ sin 9 cos 9 — — da 

d V (A3.2) 

= -\e 2f d9 + F 2 sin9 cos 9{w\ - w 2 )u\/ u 2 - I dp 

+ F 2 sin^cos^ju^ioi — w 2 )u 2 + w 2 (w 3 — w±)] d/i . 

These do in fact satisfy exactly all the d9 terms in e + and all those other terms in e + that 
have an explicit cos 9. Therefore we are left with the two remaining equations to determine the 
yet-unknown function N = N(fi,p): 

P— — = F 2 w 3 (wi — w 2 ) coshpsinhp , 

P-— = F 2 [(w 1 — w 2 ) cosh p — w\] . 
op 

We begin this part of the problem by first determining an integral for the equation involving 
ON/dp: 
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V2 N = F J dp \J (wi — w 2 ) cosh 2 p — (wi — w 3 ) 



w 3 

dp 



yj (tui - w 2 ) cosh 2 p - (wi - ty 3 ) 3 

= V^i — / <ip \l k 2 cosh 2 p — 1 Ws :F / — = 

7 V ' V™i - w s J ^/k 2 cosh 2 p-l 

a/ (wi - w 3 )E(u, k) + W3 F(u, k) 

y/Wi - W 3 

where F(u,k) and E{u,k) are the standard first and second incomplete elliptic integrals, and 
(w\ — w 2 )/(wi — w 3 ) = k 2 is the Jacobi modulus for Jacobi elliptic functions, as will be shown 
below: 

■sn _1 (z,fc) 



, IN r da r n 2 de r sn (z >y 

F(z,k)= = / = = / dw , 

(■z ./l _ ^2^2 /-sin- 1 2; pF(z,k) 

E(z,k) = da - == - = / d9\Jl-k 2 sin 2 6> = / efou dn 2 (w, fc) 
Jo VI -a 2 Jo io 



(A3.5) 



Although this is indeed the desired solution, it may be put into much more reasonable forms 
provided we now make explicit use of our solutions for the functions Wi, and also the aj, as noted 
in Eqs. (2.4) and (2.5), in terms of F = co(p + do) and the theta coefficients. However, those 
coefficients may also be expressed in terms of complete elliptic integrals: 23 

7ra 3 (/x) = 27u-^-logtf 2 (0 I ifi) = -2K(k)E(k) , 

ira 2 {n) = 2ni^- log tf 3 (0 | ifi) = -2 K(k) [E(k) - k' 2 K(k)] , 

7rai(/x) = 2iri-£- log tf 4 (0 | ifi) = -2 K(k) [E(k) - K{k)\ , ^ 3 ' 6 ^ 

Wi(fi) = ai(/j) + ^- log F(fj.) = a,i(fji) + 



dfi % n + do 

k' 2 = 1 - k 2 , n = -ir = K'(k)/K{k) , 

where the extra factor of tv appears because of our normalization, following Hancock, 23 for the 
arguments of the theta functions. Normalizations vary considerably from author to author, con- 
cerning the arguments of these functions. We will present ours at the end of this section. 

This allows us to present the two differences of functions Wi that appear in our integral in a 
different way, more easily showing the values for the integration being performed for N: 

11/1 -w 3 = -K 2 {k) , Wl -w 2 = -k 2 K 2 (k), (A3.7) 

TV TV 
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which allows for the following re-presentation of the result for N: 

F(v k) 

V^N = c (p + d ) [K(k)E(u, k) - E(k)F(u, k)} + § c Q -j^y , (A3.8) 

It is of course true that there might also be some "constant of integration," which would depend 
on 11. We show that no such constant is needed by inserting this value for N back into the 
differential equation involving its derivative with respect to n, and finding that it gives exactly 
the desired right-hand side; i.e., the pde is satisfied exactly with the value given above. 

We would like, however, to present this result also in some other formats, where the depen- 
dence on [i is made more explicit. The simplest next step is to turn it into a form involving the 
Jacobi Zeta function, Z(w, k), 23,24 and then use its relationship to the theta functions with both 
arguments non-zero. 

E(k) 

Z{w,k) = E{w,k)- 1 ^w, 
K(k)Z{F(u, k),k} = K(k)E{F(u, k), k] - E(k)F(u, k) , 
K(k)Z[2K(k)z;k] = 5^- log $±{z \ i/j) , 

K(k)Z(a, k) = rii[l, a, k] = IIi(a, k) , 
where the last line shows the relationship between the Jacobi Zeta function and the complete 
elliptic integral of the third kind, in the form originally given by Jacobi, and used by Whittaker 
and Watson. 24 When the first argument is 1 the integral is referred to as complete, and the value 
of that argument is often not shown; however, especially with Oi, the first argument is often 24 
given in terms of F(z, k) instead of z as has been done here, so that then that argument would 
be K{k). The more usual form of the integral of the third kind is the form due to Legendre, 
IL(z, a 2 , k), which was used by Olivier, 16 in his integration for the coordinates q, and q, in the 
special case when our conformal factor F is simply a constant. The Legendre form is related to 
the Jacobi form as follows, where the coefficients are various Jacobi elliptic functions, 24 and we 
also give their basic definitions, as integrals: 

dt dv 



U(z,a 2 ,k) = [ dt = f 

Jo (l-a 2 t 2 )J(l-t 2 )(l-k 2 t 2 ) Jo 



(l-a 2 t 2 )^/(l-t 2 )(l-k 2 t 2 ) Jo l-a 2 sn 2 (v,k) 

sn 2 (f, k) 
k 2 sn 2 (a, k) sn 2 (w, k) 



III {z, a, k) = k 2 sn(a, k) cn(a, k) dn(a, k) dv .„ ^ V,k) „ . - . . (A3.10) 

Jo 

ni[z, a, k] = Cn(a ' k ] ^ k) {IL[z, k 2 sn 2 (a, k),k]- z\ . 
sn(a, k) 1 
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This allows us to re- write our desired function N in several different, equivalent ways, where 
we choose the following one as most useful for our purposes: 



A/7rN(>,p) = C 



d 

tvz + + do)— log #4(2: I ip) 



(A3.ll) 

z=iF(coshp,fc)/.fs:(A;) 



where we recall here the relation of p and it to the original spherical coordinates, 6 and ip, 

m /«s •/ i sin -0 -H cos ^ cos # , .„ 

p = logtan(0/2) + , u = coshp = , (A3.12) 

sinfc' 

and, in terms of these variables, our important quantity P 2 (n,p) is given by 

P\n,p) = At- = -i [2c (n + do)K(k) dn(a, k)f . (A3.13) 
sin 9 

Completing our picture we put here the complete definitions for the theta functions, which 
are everywhere analytic functions of their first argument, z, periodic with period 1, while they 
are analytic in the upper half plane for their second argument, r. As well, they have power series 
expressions which converge very fast and are generators of the usual Jacobi elliptic functions in 
that those functions are ratios of the theta functions: 

+00 



& 4 (z\t) = l + 2^(-l) n Q n cos(2W) , 

n=l 
+00 

$ 3 (z\t) = 1 + 2^ q n2 cos(2nnz) , 

n=l 

+00 

2 (z\t) = 2 q {m+ ^ )2 cos[(2m + l)nz] , 

m=0 
+ 00 

$i(z\t) = 2 (-l) m q {rn+ ^ )2 sin[(2m + l)irz] 



(A3.14) 



m=0 



q = e inT = e- n » 



and we note again that there are various other normalizations for the arguments of these functions, 
often including the factor it into the argument, 24 so that, then, they have period n. 

A last thing to do here is to provide more detail as to how one acquires the other 3 parameters 
for the general solution 18 to the Halphen problem, via a Mobius transformation of the r (upper) 
half-plane, accompanied by appropriate transformations of the dependent functions, where in this 
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brief section we use the "overbar," as in r, to indicate the result after the transformation, rather 
than it having any relation to the complex conjugation operation: 

ar + b 



ct + d 

a +i — > a +i (T) = (cr + d) 2 a +i [T(T)] + c(ct + d) 
Wi — > Wi(j) = (cr + d) 2 Wi[r(r)] , 



> ; 



ad — be = +1 



(A3.15) 



Therefore, when we include these 3 parameters, the general solution has the following form, where 
we retain the overbars: 



~ d , n ( dr — b , 

a +i {r) = 2— log^5- 4 = + 

dr \a — cr J a — ct 

/ \ (cr + d) 2 _ _ 1 

Wi{r) = ■ h a +i {T) - c{ct + d) = o +i (r) + _ _ 

r + q t + q 



ago — b 



(A3.16) 



d - cq 

If one makes the particular, allowed choice of the 3 parameters in the transformation, of d = a = 1 
and c = = 6, then this more general form reduces to our earlier, particular form, as expected. 
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